3.1.36 \(\int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx\) [36]

3.1.36.1 Optimal result
3.1.36.2 Mathematica [A] (verified)
3.1.36.3 Rubi [A] (verified)
3.1.36.4 Maple [A] (verified)
3.1.36.5 Fricas [B] (verification not implemented)
3.1.36.6 Sympy [B] (verification not implemented)
3.1.36.7 Maxima [B] (verification not implemented)
3.1.36.8 Giac [B] (verification not implemented)
3.1.36.9 Mupad [B] (verification not implemented)

3.1.36.1 Optimal result

Integrand size = 36, antiderivative size = 156 \[ \int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx=\frac {a^2 (A+B) c^2 \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}+\frac {a^2 (3 A-8 B) c \cos ^5(e+f x)}{99 f (c-c \sin (e+f x))^7}+\frac {2 a^2 (3 A-8 B) \cos ^5(e+f x)}{693 f (c-c \sin (e+f x))^6}+\frac {2 a^2 (3 A-8 B) \cos ^5(e+f x)}{3465 c f (c-c \sin (e+f x))^5} \]

output
1/11*a^2*(A+B)*c^2*cos(f*x+e)^5/f/(c-c*sin(f*x+e))^8+1/99*a^2*(3*A-8*B)*c* 
cos(f*x+e)^5/f/(c-c*sin(f*x+e))^7+2/693*a^2*(3*A-8*B)*cos(f*x+e)^5/f/(c-c* 
sin(f*x+e))^6+2/3465*a^2*(3*A-8*B)*cos(f*x+e)^5/c/f/(c-c*sin(f*x+e))^5
 
3.1.36.2 Mathematica [A] (verified)

Time = 8.59 (sec) , antiderivative size = 285, normalized size of antiderivative = 1.83 \[ \int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx=\frac {a^2 \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right ) (1+\sin (e+f x))^2 \left (231 (27 A+28 B) \cos \left (\frac {1}{2} (e+f x)\right )-2475 (A+2 B) \cos \left (\frac {3}{2} (e+f x)\right )-2310 B \cos \left (\frac {5}{2} (e+f x)\right )-165 A \cos \left (\frac {7}{2} (e+f x)\right )+440 B \cos \left (\frac {7}{2} (e+f x)\right )+3 A \cos \left (\frac {11}{2} (e+f x)\right )-8 B \cos \left (\frac {11}{2} (e+f x)\right )+7623 A \sin \left (\frac {1}{2} (e+f x)\right )+2772 B \sin \left (\frac {1}{2} (e+f x)\right )+3465 A \sin \left (\frac {3}{2} (e+f x)\right )+2310 B \sin \left (\frac {3}{2} (e+f x)\right )-495 A \sin \left (\frac {5}{2} (e+f x)\right )-990 B \sin \left (\frac {5}{2} (e+f x)\right )+33 A \sin \left (\frac {9}{2} (e+f x)\right )-88 B \sin \left (\frac {9}{2} (e+f x)\right )\right )}{27720 c^6 f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^4 (-1+\sin (e+f x))^6} \]

input
Integrate[((a + a*Sin[e + f*x])^2*(A + B*Sin[e + f*x]))/(c - c*Sin[e + f*x 
])^6,x]
 
output
(a^2*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])*(1 + Sin[e + f*x])^2*(231*(27*A 
 + 28*B)*Cos[(e + f*x)/2] - 2475*(A + 2*B)*Cos[(3*(e + f*x))/2] - 2310*B*C 
os[(5*(e + f*x))/2] - 165*A*Cos[(7*(e + f*x))/2] + 440*B*Cos[(7*(e + f*x)) 
/2] + 3*A*Cos[(11*(e + f*x))/2] - 8*B*Cos[(11*(e + f*x))/2] + 7623*A*Sin[( 
e + f*x)/2] + 2772*B*Sin[(e + f*x)/2] + 3465*A*Sin[(3*(e + f*x))/2] + 2310 
*B*Sin[(3*(e + f*x))/2] - 495*A*Sin[(5*(e + f*x))/2] - 990*B*Sin[(5*(e + f 
*x))/2] + 33*A*Sin[(9*(e + f*x))/2] - 88*B*Sin[(9*(e + f*x))/2]))/(27720*c 
^6*f*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^4*(-1 + Sin[e + f*x])^6)
 
3.1.36.3 Rubi [A] (verified)

Time = 0.79 (sec) , antiderivative size = 149, normalized size of antiderivative = 0.96, number of steps used = 10, number of rules used = 10, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.278, Rules used = {3042, 3446, 3042, 3338, 3042, 3151, 3042, 3151, 3042, 3150}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {(a \sin (e+f x)+a)^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {(a \sin (e+f x)+a)^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6}dx\)

\(\Big \downarrow \) 3446

\(\displaystyle a^2 c^2 \int \frac {\cos ^4(e+f x) (A+B \sin (e+f x))}{(c-c \sin (e+f x))^8}dx\)

\(\Big \downarrow \) 3042

\(\displaystyle a^2 c^2 \int \frac {\cos (e+f x)^4 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^8}dx\)

\(\Big \downarrow \) 3338

\(\displaystyle a^2 c^2 \left (\frac {(3 A-8 B) \int \frac {\cos ^4(e+f x)}{(c-c \sin (e+f x))^7}dx}{11 c}+\frac {(A+B) \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}\right )\)

\(\Big \downarrow \) 3042

\(\displaystyle a^2 c^2 \left (\frac {(3 A-8 B) \int \frac {\cos (e+f x)^4}{(c-c \sin (e+f x))^7}dx}{11 c}+\frac {(A+B) \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}\right )\)

\(\Big \downarrow \) 3151

\(\displaystyle a^2 c^2 \left (\frac {(3 A-8 B) \left (\frac {2 \int \frac {\cos ^4(e+f x)}{(c-c \sin (e+f x))^6}dx}{9 c}+\frac {\cos ^5(e+f x)}{9 f (c-c \sin (e+f x))^7}\right )}{11 c}+\frac {(A+B) \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}\right )\)

\(\Big \downarrow \) 3042

\(\displaystyle a^2 c^2 \left (\frac {(3 A-8 B) \left (\frac {2 \int \frac {\cos (e+f x)^4}{(c-c \sin (e+f x))^6}dx}{9 c}+\frac {\cos ^5(e+f x)}{9 f (c-c \sin (e+f x))^7}\right )}{11 c}+\frac {(A+B) \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}\right )\)

\(\Big \downarrow \) 3151

\(\displaystyle a^2 c^2 \left (\frac {(3 A-8 B) \left (\frac {2 \left (\frac {\int \frac {\cos ^4(e+f x)}{(c-c \sin (e+f x))^5}dx}{7 c}+\frac {\cos ^5(e+f x)}{7 f (c-c \sin (e+f x))^6}\right )}{9 c}+\frac {\cos ^5(e+f x)}{9 f (c-c \sin (e+f x))^7}\right )}{11 c}+\frac {(A+B) \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}\right )\)

\(\Big \downarrow \) 3042

\(\displaystyle a^2 c^2 \left (\frac {(3 A-8 B) \left (\frac {2 \left (\frac {\int \frac {\cos (e+f x)^4}{(c-c \sin (e+f x))^5}dx}{7 c}+\frac {\cos ^5(e+f x)}{7 f (c-c \sin (e+f x))^6}\right )}{9 c}+\frac {\cos ^5(e+f x)}{9 f (c-c \sin (e+f x))^7}\right )}{11 c}+\frac {(A+B) \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}\right )\)

\(\Big \downarrow \) 3150

\(\displaystyle a^2 c^2 \left (\frac {(A+B) \cos ^5(e+f x)}{11 f (c-c \sin (e+f x))^8}+\frac {(3 A-8 B) \left (\frac {\cos ^5(e+f x)}{9 f (c-c \sin (e+f x))^7}+\frac {2 \left (\frac {\cos ^5(e+f x)}{35 c f (c-c \sin (e+f x))^5}+\frac {\cos ^5(e+f x)}{7 f (c-c \sin (e+f x))^6}\right )}{9 c}\right )}{11 c}\right )\)

input
Int[((a + a*Sin[e + f*x])^2*(A + B*Sin[e + f*x]))/(c - c*Sin[e + f*x])^6,x 
]
 
output
a^2*c^2*(((A + B)*Cos[e + f*x]^5)/(11*f*(c - c*Sin[e + f*x])^8) + ((3*A - 
8*B)*(Cos[e + f*x]^5/(9*f*(c - c*Sin[e + f*x])^7) + (2*(Cos[e + f*x]^5/(7* 
f*(c - c*Sin[e + f*x])^6) + Cos[e + f*x]^5/(35*c*f*(c - c*Sin[e + f*x])^5) 
))/(9*c)))/(11*c))
 

3.1.36.3.1 Defintions of rubi rules used

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3150
Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x 
_)])^(m_), x_Symbol] :> Simp[b*(g*Cos[e + f*x])^(p + 1)*((a + b*Sin[e + f*x 
])^m/(a*f*g*m)), x] /; FreeQ[{a, b, e, f, g, m, p}, x] && EqQ[a^2 - b^2, 0] 
 && EqQ[Simplify[m + p + 1], 0] &&  !ILtQ[p, 0]
 

rule 3151
Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x 
_)])^(m_), x_Symbol] :> Simp[b*(g*Cos[e + f*x])^(p + 1)*((a + b*Sin[e + f*x 
])^m/(a*f*g*Simplify[2*m + p + 1])), x] + Simp[Simplify[m + p + 1]/(a*Simpl 
ify[2*m + p + 1])   Int[(g*Cos[e + f*x])^p*(a + b*Sin[e + f*x])^(m + 1), x] 
, x] /; FreeQ[{a, b, e, f, g, m, p}, x] && EqQ[a^2 - b^2, 0] && ILtQ[Simpli 
fy[m + p + 1], 0] && NeQ[2*m + p + 1, 0] &&  !IGtQ[m, 0]
 

rule 3338
Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x 
_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Simp[(b*c - 
 a*d)*(g*Cos[e + f*x])^(p + 1)*((a + b*Sin[e + f*x])^m/(a*f*g*(2*m + p + 1) 
)), x] + Simp[(a*d*m + b*c*(m + p + 1))/(a*b*(2*m + p + 1))   Int[(g*Cos[e 
+ f*x])^p*(a + b*Sin[e + f*x])^(m + 1), x], x] /; FreeQ[{a, b, c, d, e, f, 
g, m, p}, x] && EqQ[a^2 - b^2, 0] && (LtQ[m, -1] || ILtQ[Simplify[m + p], 0 
]) && NeQ[2*m + p + 1, 0]
 

rule 3446
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) + 
 (f_.)*(x_)])*((c_) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Si 
mp[a^m*c^m   Int[Cos[e + f*x]^(2*m)*(c + d*Sin[e + f*x])^(n - m)*(A + B*Sin 
[e + f*x]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B, n}, x] && EqQ[b*c + a* 
d, 0] && EqQ[a^2 - b^2, 0] && IntegerQ[m] &&  !(IntegerQ[n] && ((LtQ[m, 0] 
&& GtQ[n, 0]) || LtQ[0, n, m] || LtQ[m, n, 0]))
 
3.1.36.4 Maple [A] (verified)

Time = 1.49 (sec) , antiderivative size = 219, normalized size of antiderivative = 1.40

method result size
parallelrisch \(-\frac {2 a^{2} \left (A \left (\tan ^{10}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\left (-3 A +B \right ) \left (\tan ^{9}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\left (12 A -\frac {B}{3}\right ) \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+2 \left (-10 A +\frac {7 B}{3}\right ) \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )+\frac {2 \left (81 A -\frac {13 B}{3}\right ) \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{5}+\frac {2 \left (-71 A +16 B \right ) \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{5}+\frac {4 \left (41 A -2 B \right ) \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{7}+\frac {2 \left (-34 A +9 B \right ) \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{7}+\frac {\left (89 A +\frac {2 B}{3}\right ) \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{21}+\frac {\left (-47 A +\frac {61 B}{3}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{105}+\frac {152 A}{1155}-\frac {61 B}{3465}\right )}{f \,c^{6} \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{11}}\) \(219\)
derivativedivides \(\frac {2 a^{2} \left (-\frac {1752 A +1208 B}{6 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{6}}-\frac {A}{\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1}-\frac {90 A +26 B}{3 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{3}}-\frac {640 A +640 B}{10 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{10}}-\frac {2304 A +2048 B}{8 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{8}}-\frac {128 A +128 B}{11 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{11}}-\frac {1536 A +1472 B}{9 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{9}}-\frac {14 A +2 B}{2 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{2}}-\frac {932 A +528 B}{5 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{5}}-\frac {2376 A +1896 B}{7 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{7}}-\frac {352 A +152 B}{4 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{4}}\right )}{f \,c^{6}}\) \(249\)
default \(\frac {2 a^{2} \left (-\frac {1752 A +1208 B}{6 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{6}}-\frac {A}{\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1}-\frac {90 A +26 B}{3 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{3}}-\frac {640 A +640 B}{10 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{10}}-\frac {2304 A +2048 B}{8 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{8}}-\frac {128 A +128 B}{11 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{11}}-\frac {1536 A +1472 B}{9 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{9}}-\frac {14 A +2 B}{2 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{2}}-\frac {932 A +528 B}{5 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{5}}-\frac {2376 A +1896 B}{7 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{7}}-\frac {352 A +152 B}{4 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{4}}\right )}{f \,c^{6}}\) \(249\)
risch \(\frac {-\frac {4 A \,a^{2}}{1155}-\frac {16 i B \,a^{2} {\mathrm e}^{5 i \left (f x +e \right )}}{5}-\frac {44 i A \,a^{2} {\mathrm e}^{5 i \left (f x +e \right )}}{5}+\frac {32 B \,a^{2}}{3465}+\frac {32 i B \,a^{2} {\mathrm e}^{i \left (f x +e \right )}}{315}-\frac {4 i A \,a^{2} {\mathrm e}^{i \left (f x +e \right )}}{105}+\frac {8 i B \,a^{2} {\mathrm e}^{7 i \left (f x +e \right )}}{3}-\frac {32 B \,a^{2} {\mathrm e}^{2 i \left (f x +e \right )}}{63}+4 i A \,a^{2} {\mathrm e}^{7 i \left (f x +e \right )}+\frac {8 i B \,a^{2} {\mathrm e}^{3 i \left (f x +e \right )}}{7}+\frac {4 i A \,a^{2} {\mathrm e}^{3 i \left (f x +e \right )}}{7}+\frac {8 B \,a^{2} {\mathrm e}^{8 i \left (f x +e \right )}}{3}-\frac {36 A \,a^{2} {\mathrm e}^{6 i \left (f x +e \right )}}{5}-\frac {112 B \,a^{2} {\mathrm e}^{6 i \left (f x +e \right )}}{15}+\frac {20 A \,a^{2} {\mathrm e}^{4 i \left (f x +e \right )}}{7}+\frac {40 B \,a^{2} {\mathrm e}^{4 i \left (f x +e \right )}}{7}+\frac {4 A \,a^{2} {\mathrm e}^{2 i \left (f x +e \right )}}{21}}{\left ({\mathrm e}^{i \left (f x +e \right )}-i\right )^{11} f \,c^{6}}\) \(269\)
norman \(\frac {-\frac {912 A \,a^{2}-122 B \,a^{2}}{3465 c f}-\frac {2 A \,a^{2} \left (\tan ^{16}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c f}+\frac {2 \left (3 A \,a^{2}-B \,a^{2}\right ) \left (\tan ^{15}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c f}-\frac {2 \left (45 A \,a^{2}-B \,a^{2}\right ) \left (\tan ^{14}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 c f}+\frac {2 \left (87 A \,a^{2}-23 B \,a^{2}\right ) \left (\tan ^{13}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 c f}+\frac {\left (282 A \,a^{2}-122 B \,a^{2}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{315 c f}+\frac {2 \left (487 A \,a^{2}-117 B \,a^{2}\right ) \left (\tan ^{11}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{5 c f}-\frac {2 \left (1071 A \,a^{2}-41 B \,a^{2}\right ) \left (\tan ^{12}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{15 c f}+\frac {2 \left (1161 A \,a^{2}-331 B \,a^{2}\right ) \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{105 c f}-\frac {2 \left (5517 A \,a^{2}-257 B \,a^{2}\right ) \left (\tan ^{10}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{35 c f}+\frac {2 \left (5527 A \,a^{2}-1287 B \,a^{2}\right ) \left (\tan ^{9}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{35 c f}+\frac {2 \left (6183 A \,a^{2}-1543 B \,a^{2}\right ) \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{105 c f}-\frac {2 \left (16053 A \,a^{2}-73 B \,a^{2}\right ) \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3465 c f}-\frac {2 \left (42201 A \,a^{2}-1271 B \,a^{2}\right ) \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{1155 c f}+\frac {2 \left (42459 A \,a^{2}-10009 B \,a^{2}\right ) \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{315 c f}-\frac {2 \left (57873 A \,a^{2}-2813 B \,a^{2}\right ) \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{315 c f}-\frac {2 \left (400317 A \,a^{2}-17617 B \,a^{2}\right ) \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3465 c f}}{\left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{3} c^{5} \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{11}}\) \(556\)

input
int((a+a*sin(f*x+e))^2*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^6,x,method=_RETUR 
NVERBOSE)
 
output
-2*a^2*(A*tan(1/2*f*x+1/2*e)^10+(-3*A+B)*tan(1/2*f*x+1/2*e)^9+(12*A-1/3*B) 
*tan(1/2*f*x+1/2*e)^8+2*(-10*A+7/3*B)*tan(1/2*f*x+1/2*e)^7+2/5*(81*A-13/3* 
B)*tan(1/2*f*x+1/2*e)^6+2/5*(-71*A+16*B)*tan(1/2*f*x+1/2*e)^5+4/7*(41*A-2* 
B)*tan(1/2*f*x+1/2*e)^4+2/7*(-34*A+9*B)*tan(1/2*f*x+1/2*e)^3+1/21*(89*A+2/ 
3*B)*tan(1/2*f*x+1/2*e)^2+1/105*(-47*A+61/3*B)*tan(1/2*f*x+1/2*e)+152/1155 
*A-61/3465*B)/f/c^6/(tan(1/2*f*x+1/2*e)-1)^11
 
3.1.36.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 407 vs. \(2 (152) = 304\).

Time = 0.25 (sec) , antiderivative size = 407, normalized size of antiderivative = 2.61 \[ \int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx=-\frac {2 \, {\left (3 \, A - 8 \, B\right )} a^{2} \cos \left (f x + e\right )^{6} + 12 \, {\left (3 \, A - 8 \, B\right )} a^{2} \cos \left (f x + e\right )^{5} - 25 \, {\left (3 \, A - 8 \, B\right )} a^{2} \cos \left (f x + e\right )^{4} - 35 \, {\left (6 \, A + 17 \, B\right )} a^{2} \cos \left (f x + e\right )^{3} - 35 \, {\left (21 \, A + 43 \, B\right )} a^{2} \cos \left (f x + e\right )^{2} + 630 \, {\left (A + B\right )} a^{2} \cos \left (f x + e\right ) + 1260 \, {\left (A + B\right )} a^{2} - {\left (2 \, {\left (3 \, A - 8 \, B\right )} a^{2} \cos \left (f x + e\right )^{5} - 10 \, {\left (3 \, A - 8 \, B\right )} a^{2} \cos \left (f x + e\right )^{4} - 35 \, {\left (3 \, A - 8 \, B\right )} a^{2} \cos \left (f x + e\right )^{3} + 35 \, {\left (3 \, A + 25 \, B\right )} a^{2} \cos \left (f x + e\right )^{2} - 630 \, {\left (A + B\right )} a^{2} \cos \left (f x + e\right ) - 1260 \, {\left (A + B\right )} a^{2}\right )} \sin \left (f x + e\right )}{3465 \, {\left (c^{6} f \cos \left (f x + e\right )^{6} - 5 \, c^{6} f \cos \left (f x + e\right )^{5} - 18 \, c^{6} f \cos \left (f x + e\right )^{4} + 20 \, c^{6} f \cos \left (f x + e\right )^{3} + 48 \, c^{6} f \cos \left (f x + e\right )^{2} - 16 \, c^{6} f \cos \left (f x + e\right ) - 32 \, c^{6} f + {\left (c^{6} f \cos \left (f x + e\right )^{5} + 6 \, c^{6} f \cos \left (f x + e\right )^{4} - 12 \, c^{6} f \cos \left (f x + e\right )^{3} - 32 \, c^{6} f \cos \left (f x + e\right )^{2} + 16 \, c^{6} f \cos \left (f x + e\right ) + 32 \, c^{6} f\right )} \sin \left (f x + e\right )\right )}} \]

input
integrate((a+a*sin(f*x+e))^2*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^6,x, algori 
thm="fricas")
 
output
-1/3465*(2*(3*A - 8*B)*a^2*cos(f*x + e)^6 + 12*(3*A - 8*B)*a^2*cos(f*x + e 
)^5 - 25*(3*A - 8*B)*a^2*cos(f*x + e)^4 - 35*(6*A + 17*B)*a^2*cos(f*x + e) 
^3 - 35*(21*A + 43*B)*a^2*cos(f*x + e)^2 + 630*(A + B)*a^2*cos(f*x + e) + 
1260*(A + B)*a^2 - (2*(3*A - 8*B)*a^2*cos(f*x + e)^5 - 10*(3*A - 8*B)*a^2* 
cos(f*x + e)^4 - 35*(3*A - 8*B)*a^2*cos(f*x + e)^3 + 35*(3*A + 25*B)*a^2*c 
os(f*x + e)^2 - 630*(A + B)*a^2*cos(f*x + e) - 1260*(A + B)*a^2)*sin(f*x + 
 e))/(c^6*f*cos(f*x + e)^6 - 5*c^6*f*cos(f*x + e)^5 - 18*c^6*f*cos(f*x + e 
)^4 + 20*c^6*f*cos(f*x + e)^3 + 48*c^6*f*cos(f*x + e)^2 - 16*c^6*f*cos(f*x 
 + e) - 32*c^6*f + (c^6*f*cos(f*x + e)^5 + 6*c^6*f*cos(f*x + e)^4 - 12*c^6 
*f*cos(f*x + e)^3 - 32*c^6*f*cos(f*x + e)^2 + 16*c^6*f*cos(f*x + e) + 32*c 
^6*f)*sin(f*x + e))
 
3.1.36.6 Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4816 vs. \(2 (141) = 282\).

Time = 44.16 (sec) , antiderivative size = 4816, normalized size of antiderivative = 30.87 \[ \int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx=\text {Too large to display} \]

input
integrate((a+a*sin(f*x+e))**2*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))**6,x)
 
output
Piecewise((-6930*A*a**2*tan(e/2 + f*x/2)**10/(3465*c**6*f*tan(e/2 + f*x/2) 
**11 - 38115*c**6*f*tan(e/2 + f*x/2)**10 + 190575*c**6*f*tan(e/2 + f*x/2)* 
*9 - 571725*c**6*f*tan(e/2 + f*x/2)**8 + 1143450*c**6*f*tan(e/2 + f*x/2)** 
7 - 1600830*c**6*f*tan(e/2 + f*x/2)**6 + 1600830*c**6*f*tan(e/2 + f*x/2)** 
5 - 1143450*c**6*f*tan(e/2 + f*x/2)**4 + 571725*c**6*f*tan(e/2 + f*x/2)**3 
 - 190575*c**6*f*tan(e/2 + f*x/2)**2 + 38115*c**6*f*tan(e/2 + f*x/2) - 346 
5*c**6*f) + 20790*A*a**2*tan(e/2 + f*x/2)**9/(3465*c**6*f*tan(e/2 + f*x/2) 
**11 - 38115*c**6*f*tan(e/2 + f*x/2)**10 + 190575*c**6*f*tan(e/2 + f*x/2)* 
*9 - 571725*c**6*f*tan(e/2 + f*x/2)**8 + 1143450*c**6*f*tan(e/2 + f*x/2)** 
7 - 1600830*c**6*f*tan(e/2 + f*x/2)**6 + 1600830*c**6*f*tan(e/2 + f*x/2)** 
5 - 1143450*c**6*f*tan(e/2 + f*x/2)**4 + 571725*c**6*f*tan(e/2 + f*x/2)**3 
 - 190575*c**6*f*tan(e/2 + f*x/2)**2 + 38115*c**6*f*tan(e/2 + f*x/2) - 346 
5*c**6*f) - 83160*A*a**2*tan(e/2 + f*x/2)**8/(3465*c**6*f*tan(e/2 + f*x/2) 
**11 - 38115*c**6*f*tan(e/2 + f*x/2)**10 + 190575*c**6*f*tan(e/2 + f*x/2)* 
*9 - 571725*c**6*f*tan(e/2 + f*x/2)**8 + 1143450*c**6*f*tan(e/2 + f*x/2)** 
7 - 1600830*c**6*f*tan(e/2 + f*x/2)**6 + 1600830*c**6*f*tan(e/2 + f*x/2)** 
5 - 1143450*c**6*f*tan(e/2 + f*x/2)**4 + 571725*c**6*f*tan(e/2 + f*x/2)**3 
 - 190575*c**6*f*tan(e/2 + f*x/2)**2 + 38115*c**6*f*tan(e/2 + f*x/2) - 346 
5*c**6*f) + 138600*A*a**2*tan(e/2 + f*x/2)**7/(3465*c**6*f*tan(e/2 + f*x/2 
)**11 - 38115*c**6*f*tan(e/2 + f*x/2)**10 + 190575*c**6*f*tan(e/2 + f*x...
 
3.1.36.7 Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 2604 vs. \(2 (152) = 304\).

Time = 0.31 (sec) , antiderivative size = 2604, normalized size of antiderivative = 16.69 \[ \int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx=\text {Too large to display} \]

input
integrate((a+a*sin(f*x+e))^2*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^6,x, algori 
thm="maxima")
 
output
-2/3465*(5*A*a^2*(913*sin(f*x + e)/(cos(f*x + e) + 1) - 4565*sin(f*x + e)^ 
2/(cos(f*x + e) + 1)^2 + 12540*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 25080 
*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 33726*sin(f*x + e)^5/(cos(f*x + e) 
+ 1)^5 - 33726*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 + 23100*sin(f*x + e)^7/ 
(cos(f*x + e) + 1)^7 - 11550*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 + 3465*si 
n(f*x + e)^9/(cos(f*x + e) + 1)^9 - 693*sin(f*x + e)^10/(cos(f*x + e) + 1) 
^10 - 146)/(c^6 - 11*c^6*sin(f*x + e)/(cos(f*x + e) + 1) + 55*c^6*sin(f*x 
+ e)^2/(cos(f*x + e) + 1)^2 - 165*c^6*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 
+ 330*c^6*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - 462*c^6*sin(f*x + e)^5/(co 
s(f*x + e) + 1)^5 + 462*c^6*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - 330*c^6* 
sin(f*x + e)^7/(cos(f*x + e) + 1)^7 + 165*c^6*sin(f*x + e)^8/(cos(f*x + e) 
 + 1)^8 - 55*c^6*sin(f*x + e)^9/(cos(f*x + e) + 1)^9 + 11*c^6*sin(f*x + e) 
^10/(cos(f*x + e) + 1)^10 - c^6*sin(f*x + e)^11/(cos(f*x + e) + 1)^11) - 6 
*A*a^2*(671*sin(f*x + e)/(cos(f*x + e) + 1) - 2200*sin(f*x + e)^2/(cos(f*x 
 + e) + 1)^2 + 6600*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 10890*sin(f*x + 
e)^4/(cos(f*x + e) + 1)^4 + 15246*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 - 12 
936*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 + 9240*sin(f*x + e)^7/(cos(f*x + e 
) + 1)^7 - 3465*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 + 1155*sin(f*x + e)^9/ 
(cos(f*x + e) + 1)^9 - 61)/(c^6 - 11*c^6*sin(f*x + e)/(cos(f*x + e) + 1) + 
 55*c^6*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 165*c^6*sin(f*x + e)^3/(c...
 
3.1.36.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 353 vs. \(2 (152) = 304\).

Time = 0.36 (sec) , antiderivative size = 353, normalized size of antiderivative = 2.26 \[ \int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx=-\frac {2 \, {\left (3465 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{10} - 10395 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{9} + 3465 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{9} + 41580 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{8} - 1155 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{8} - 69300 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{7} + 16170 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{7} + 112266 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{6} - 6006 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{6} - 98406 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{5} + 22176 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{5} + 81180 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{4} - 3960 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{4} - 33660 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} + 8910 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} + 14685 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 110 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - 1551 \, A a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 671 \, B a^{2} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 456 \, A a^{2} - 61 \, B a^{2}\right )}}{3465 \, c^{6} f {\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 1\right )}^{11}} \]

input
integrate((a+a*sin(f*x+e))^2*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^6,x, algori 
thm="giac")
 
output
-2/3465*(3465*A*a^2*tan(1/2*f*x + 1/2*e)^10 - 10395*A*a^2*tan(1/2*f*x + 1/ 
2*e)^9 + 3465*B*a^2*tan(1/2*f*x + 1/2*e)^9 + 41580*A*a^2*tan(1/2*f*x + 1/2 
*e)^8 - 1155*B*a^2*tan(1/2*f*x + 1/2*e)^8 - 69300*A*a^2*tan(1/2*f*x + 1/2* 
e)^7 + 16170*B*a^2*tan(1/2*f*x + 1/2*e)^7 + 112266*A*a^2*tan(1/2*f*x + 1/2 
*e)^6 - 6006*B*a^2*tan(1/2*f*x + 1/2*e)^6 - 98406*A*a^2*tan(1/2*f*x + 1/2* 
e)^5 + 22176*B*a^2*tan(1/2*f*x + 1/2*e)^5 + 81180*A*a^2*tan(1/2*f*x + 1/2* 
e)^4 - 3960*B*a^2*tan(1/2*f*x + 1/2*e)^4 - 33660*A*a^2*tan(1/2*f*x + 1/2*e 
)^3 + 8910*B*a^2*tan(1/2*f*x + 1/2*e)^3 + 14685*A*a^2*tan(1/2*f*x + 1/2*e) 
^2 + 110*B*a^2*tan(1/2*f*x + 1/2*e)^2 - 1551*A*a^2*tan(1/2*f*x + 1/2*e) + 
671*B*a^2*tan(1/2*f*x + 1/2*e) + 456*A*a^2 - 61*B*a^2)/(c^6*f*(tan(1/2*f*x 
 + 1/2*e) - 1)^11)
 
3.1.36.9 Mupad [B] (verification not implemented)

Time = 15.72 (sec) , antiderivative size = 423, normalized size of antiderivative = 2.71 \[ \int \frac {(a+a \sin (e+f x))^2 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx=\frac {2\,\cos \left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (\frac {38163\,A\,a^2}{8}-\frac {1283\,B\,a^2}{8}-\frac {11931\,A\,a^2\,\cos \left (2\,e+2\,f\,x\right )}{4}+\frac {9609\,A\,a^2\,\cos \left (3\,e+3\,f\,x\right )}{16}+\frac {1383\,A\,a^2\,\cos \left (4\,e+4\,f\,x\right )}{8}-\frac {225\,A\,a^2\,\cos \left (5\,e+5\,f\,x\right )}{16}+\frac {631\,B\,a^2\,\cos \left (2\,e+2\,f\,x\right )}{4}-\frac {1583\,B\,a^2\,\cos \left (3\,e+3\,f\,x\right )}{32}-\frac {223\,B\,a^2\,\cos \left (4\,e+4\,f\,x\right )}{8}+\frac {45\,B\,a^2\,\cos \left (5\,e+5\,f\,x\right )}{32}+1386\,A\,a^2\,\sin \left (2\,e+2\,f\,x\right )+\frac {14949\,A\,a^2\,\sin \left (3\,e+3\,f\,x\right )}{16}-\frac {561\,A\,a^2\,\sin \left (4\,e+4\,f\,x\right )}{4}-\frac {231\,A\,a^2\,\sin \left (5\,e+5\,f\,x\right )}{16}-\frac {3003\,B\,a^2\,\sin \left (2\,e+2\,f\,x\right )}{8}-\frac {4653\,B\,a^2\,\sin \left (3\,e+3\,f\,x\right )}{32}+\frac {209\,B\,a^2\,\sin \left (4\,e+4\,f\,x\right )}{16}+\frac {77\,B\,a^2\,\sin \left (5\,e+5\,f\,x\right )}{32}-2091\,A\,a^2\,\cos \left (e+f\,x\right )+\frac {281\,B\,a^2\,\cos \left (e+f\,x\right )}{16}-\frac {22869\,A\,a^2\,\sin \left (e+f\,x\right )}{4}+\frac {23331\,B\,a^2\,\sin \left (e+f\,x\right )}{16}\right )}{3465\,c^6\,f\,\left (\frac {231\,\sqrt {2}\,\cos \left (\frac {e}{2}+\frac {\pi }{4}+\frac {f\,x}{2}\right )}{16}-\frac {165\,\sqrt {2}\,\cos \left (\frac {3\,e}{2}-\frac {\pi }{4}+\frac {3\,f\,x}{2}\right )}{16}-\frac {165\,\sqrt {2}\,\cos \left (\frac {5\,e}{2}+\frac {\pi }{4}+\frac {5\,f\,x}{2}\right )}{32}+\frac {55\,\sqrt {2}\,\cos \left (\frac {7\,e}{2}-\frac {\pi }{4}+\frac {7\,f\,x}{2}\right )}{32}+\frac {11\,\sqrt {2}\,\cos \left (\frac {9\,e}{2}+\frac {\pi }{4}+\frac {9\,f\,x}{2}\right )}{32}-\frac {\sqrt {2}\,\cos \left (\frac {11\,e}{2}-\frac {\pi }{4}+\frac {11\,f\,x}{2}\right )}{32}\right )} \]

input
int(((A + B*sin(e + f*x))*(a + a*sin(e + f*x))^2)/(c - c*sin(e + f*x))^6,x 
)
 
output
(2*cos(e/2 + (f*x)/2)*((38163*A*a^2)/8 - (1283*B*a^2)/8 - (11931*A*a^2*cos 
(2*e + 2*f*x))/4 + (9609*A*a^2*cos(3*e + 3*f*x))/16 + (1383*A*a^2*cos(4*e 
+ 4*f*x))/8 - (225*A*a^2*cos(5*e + 5*f*x))/16 + (631*B*a^2*cos(2*e + 2*f*x 
))/4 - (1583*B*a^2*cos(3*e + 3*f*x))/32 - (223*B*a^2*cos(4*e + 4*f*x))/8 + 
 (45*B*a^2*cos(5*e + 5*f*x))/32 + 1386*A*a^2*sin(2*e + 2*f*x) + (14949*A*a 
^2*sin(3*e + 3*f*x))/16 - (561*A*a^2*sin(4*e + 4*f*x))/4 - (231*A*a^2*sin( 
5*e + 5*f*x))/16 - (3003*B*a^2*sin(2*e + 2*f*x))/8 - (4653*B*a^2*sin(3*e + 
 3*f*x))/32 + (209*B*a^2*sin(4*e + 4*f*x))/16 + (77*B*a^2*sin(5*e + 5*f*x) 
)/32 - 2091*A*a^2*cos(e + f*x) + (281*B*a^2*cos(e + f*x))/16 - (22869*A*a^ 
2*sin(e + f*x))/4 + (23331*B*a^2*sin(e + f*x))/16))/(3465*c^6*f*((231*2^(1 
/2)*cos(e/2 + pi/4 + (f*x)/2))/16 - (165*2^(1/2)*cos((3*e)/2 - pi/4 + (3*f 
*x)/2))/16 - (165*2^(1/2)*cos((5*e)/2 + pi/4 + (5*f*x)/2))/32 + (55*2^(1/2 
)*cos((7*e)/2 - pi/4 + (7*f*x)/2))/32 + (11*2^(1/2)*cos((9*e)/2 + pi/4 + ( 
9*f*x)/2))/32 - (2^(1/2)*cos((11*e)/2 - pi/4 + (11*f*x)/2))/32))